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Abstract 

A model of laminated wave turbulence is presented. This model consists of two 
co-existing layers - one with continuous waves' spectra, covered by KAM theory and 
Kolmogorov-like power spectra, and one with discrete waves' spectra, covered by discrete 
classes of waves and Clipping method. Some known laboratory experiments and numeri- 
cal simulations are explained in the frame of this model. 

PACS: 47.10.-g, 47.27.De, 47.27.T 

1. Continuous wave spectra 

In ^2] Kolmogorov presented energy spectrum of turbulence describing the distribu- 
tion of the energy among turbulence vortices as function of vortex size and thus founded 
the field of mathematical analysis of turbulence. Kolmogorov regarded some inertial range 
of wave numbers, between viscosity and dissipation, and suggested that at this range, tur- 
bulence is (1) locally homogeneous (no dependence on position) and (2) locally isotropic 
(no dependence on direction) which can be summarized as follows: probability distribution 
for the relative velocities of two particles in the fluid only depends on the distance between 
particles. Using these suggestions and dimensional analysis, Kolmogorov deduced that en- 
ergy distribution, called now Kolmogorov 's spectrum, is proportional to A;~ 5 / 3 for wave 
numbers k. Results of numerical simulations and real experiments carried out to prove 
this theory are somewhat contradictious. On the one hand, probably the most spectacular 
example of the validity of Kolmogorov 's spectra is provided in [3] where measurements in 
tidal currents near Seymour Narrows north of Campbell River on Vancouver Island were 
described and —5/3 spectra appeared at the range of 10 4 (energy dissipation at a scale of 
millimeters and energy input - at 100 meters). On the other hand, Kolmogorov 's spectra 
have been obtained under the assumptions opposite to Kolmogorov 's Hj so that exponent 
—5/3 corresponds to both direct and inverse cascades. 

With a hope to diminish established unclearness of Kolmogorov 's theory in a more 
simple setting, theory of wave (or weak) turbulence (WT) in the systems with continuous 
wave spectra has been developed. The problem is regarded in the very general form 

LU>) = sN(ij) (1) 

where L and N are linear and nonlinear operators consequently, 0<e<<lisa parameter 
of nonlinearity and linear part possesses wave-like solutions of the form 

ijj(x) = A(k, x) exp i(kx — uj(k)t). (2) 
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Choice of e is defined by specifics of the physical wave system under the study. For 
instance, for spherical planetary waves e is usually chosen as the ratio of the particle 
velocity to the phase velocity while for water waves usually e = a\k\ , where a is amplitude 
of a wave and e characterizes in this case the steepness of the waves. For small enough e, 
solutions of Eq.Q are described at the slow-time scales T\ = t/e, T 2 = t/e 2 , T3 = i/e 3 , 
etc. by resonantly interacting waves only, i.e. by the waves with wave-vectors satisfying 
to resonant conditions (for nonlinearity of order n): 

w(fci) ± io(k 2 ) ± ... ± oj(k n+1 ) = 0, 
ki ± k 2 ± ... ± k n +i = 

so that quadratic nonlinearity corresponds to 3- waves interactions, cubic - to 4- waves 
interactions, etc. Obviously, in every physical problem the resonances have some nonzero 
width, i.e. Eq.((3Jl) takes form 

io(ki) ± oj(k 2 ) ± ... ± u(k n+ i) = A (4) 

with a small but nonzero discrepancy < A << 1. Taking this into account, nonlinear 
part of Q can be rewritten as 

v Vi8{k ± k 2 ± ... ± ki) 
Lj(k%) ± oj{k 2 ) ± ... ± w(fcj) 

where <5 is a delta-function and each term of this sum with nonzero vertex coefficient 
Vi corresponds to a specific slow-time scale of wave interactions. The representation 
(jSJ is used then for construction of a wave kinetic equation, with corresponding vertex 
coefficients and delta-functions in the under-integral expression. 

From mathematical point of view, it is important to establish the finiteness of non- 
linearity given by © in the case when A — ► because representation © becomes mean- 
ingless if A = 0. This problem - so-called "problem of small denominators" - was solved 
by KAM theory (|14j.[T].|16j) in the following way. For small enough A and sufficiently 
irrational dispersion function u, these wave systems contain an infinite set of invariant tori 
which carry quasi-periodic motions which in phase space are confined to the tori. Main 
result of KAM theory is therefore a decomposition of action into disjoint invariant sets, 
and though it contradicts ergodicity but not very substantially as the size of the system 
tends to infinity 2 . In particularly, random phase approximation can be assumed and 
kinetic equations and Kohnogorov-like power spectra /c 7 , 7 < 0, give then appropriate 
description of these wave systems at the corresponding time scales. It means that some 
special set of points in spectral space, corresponding to the A-vicinity of exact resonances, 
has been excluded from consideration in order to obtain wave kinetic equation. We give 
detailed description of this subset of spectral space in the section 4 of this paper. 

2. Discrete wave spectra 

There exist a lot of wave phenomena which are due to discreteness of the wave spectra 
(corresponds to zero- or periodic boundary conditions) and can not be explained in terms 
of kinetic equations and power energy spectra. To describe these phenomena, WT in the 
systems with discrete spectra has been developed ^0])CB- ^ turned out that discrete 
systems possess some qualitatively new properties: (a) all resonantly interacting waves 
are divided into disjoint classes, there is no energy flow between different these classes; (b) 
major part of the waves do not interact; (c) all interactions of a specific wave are confined 
to some finite domain; (d) number of interacting waves depends on the form of boundary 
conditions, for the great number of boundary conditions interactions are not possible; (e) 
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all properties (a)-(d) keep true for approximate interactions, i.e. for some small enough 
discrepancy < A << 1. This fact gave a rise to Clipping method 9 j which allows "to 
clip out" all non-interacting waves from the whole spectra and study only those which do 
interact, exactly or approximately. Approximate resonances are understood on a discrete 
lattice, i.e. wave vectors of approximately interacting waves are also integers. The ener- 
getic behavior of these systems is described then independently (at each slow-time scale 
Tj) by a few small systems of ordinary differential equations (SODE) on slowly-changing 
amplitudes of resonantly interacting waves, i.e. amplitude of linear wave (J2J) is a function 
of some Tj depending on the form on nonlinearity N. For instance, for 3-waves interactions 
major part of these SODE consist of three equations on three (real-valued) amplitudes 
and can be solved explicitly in terms of elliptic functions on T\. To compare with WT 
of the systems with continuous spectra, SODE are to be used instead of kinetic equation 
and their coupling coefficients - instead of power-law spectra. Notice that though kinetic 
description does not apply for discrete systems, some of these results are in a sense similar 
to those of KAM theory, for instance, Theorem on the partition [Jj can be regarded as an 
analog of KAM-Theorem for discrete systems. 

3. Transition from discrete to continuous spectra 

Now, the standard qualitative model of the WT can be presented as follows: short 
waves are described by Kolmogorov's energy spectra and kinetic equations, long waves 
are described by Clipping method and dynamic equations, and somewhere in between a 
"transition" interval exists that has its own specifics and should be described separately. 

We would like to demonstrate some contradictiveness of this qualitative model and 
begin with two remarks. (1) WT of discrete waves systems has been developed for ar- 
bitrary wave numbers which means that transition from finite to infinite domain can be 
constructed not only in some finite "transition" interval but at the whole infinite range 
of wave numbers. (2) Transition from discrete to continuous spectrum is often regarded 
in somewhat over-simplified way: if say, real- valued wave vectors k = (k x ,k y ) have dis- 
persion function uj(k x ,k y ) with k x ,k y £ R, then the same function of integer variables, 
u(m, n) with m, n £ Z, describes corresponding discrete waves. In general, it is not true. 
We demonstrate it taking barotropic vorticity equation (BVE), also known as Obukhov- 
Charney-Hasegawa-Mima equation, as our main example motivated by its wide applica- 
bility for describing a great number of physically important phenomena in astrophysics, 
geophysics and plasma physics. 

BVE on a sphere has form 

^|* + 2^ + J W ,A« = (0) 

with linear waves of the form 

^sphere = AP™ (sin (f>) exp i[m\ + ^sphere t]. (7) 

Here ip is the stream-function; variables t, (ft and A physically mean the time, the latitude 
(— 7r/2 < eft < it/2) and the longitude (0 < A < 2ir) respectively; P™(x) is the associated 
Legendre function of degree n and order m. The same equation taken on infinite /3-plane 
has linear waves of the form 

l^plane = A exp i(k x X + kylj + iO p i ane t), 

which means that 

^sphere = m/[n(n + 1)] and UJplane = K/ (1 + kl + ky), 
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here constant multipliers are omitted because they disappear due to homogeneous form 
of Eq.(J31l). It is easy to see that no wave vectors k = (m,n) : m,n G Z satisfy Eq.^l) 
with oJsphere an d with Lu p i ane simultaneously. It means that discrete waves do not have 
images on infinite plane when such a "naive" transition is regarded. 

More intrinsic construction of the transition from spherical to plane planetary waves 
IHj can be derived in following way. Regarding m ~ n >> 1 and using asymptotic 
approximation for Legendre functions, one can "convert" (not always but in a bounded 
latitudinal belt with the width ~ n _1 ) one spherical wave into a linear combination of 
two plane waves 

Aexpi(k((po) x x ± k((p ) y y + UJ p l ane t), 

where local wave numbers k((fo) x , k(ipo) y £ R are functions of the initial spherical wave 
number m,n and of the so-called interaction latitude c^o : 

2 m 2 (n 2 . + n| — n 2 ) + m\{p\ + n 2 — n 2 ,) + m^(n 2 + n 2 , — n^)) 

n\ n 2 , + nfnl + n^n 2 . — (n\ + n\ + n|)/4 

If interaction latitude exists, < cos 2 (po < 1, plane images of spherical waves interact as 
in classical /3-plane approximation. In particularly, this means that (1) transition from a 
spherical domain to an infinite plane is transition to a one-parametric family of infinite 
planes, and (2) such a transition is not always possible. A very important fact is that 
plane wave system keeps memory about spherical interactions: coupling coefficient of the 
plane images of spherical waves is ~ k 3 / 2 and ~ A; 7 / 6 otherwise and k = \k\. 

The same reasoning allows to construct a transition from a square domain (dispersion 
function being then u! sqU are = l/Vm 2 + n 2 ) to the infinite /?-plane where difference in 
magnitudes of coupling coefficients is even more substantial: ~ k 2 for plane images of the 
waves from square domain and ~ k otherwise 0. These results hold for discrete approx- 
imate interactions in following way: long-wave part of spectrum is dominated by a few 
resonantly interacting waves with huge amplitudes while short-wave part of the spectrum 
consists of many approximately interacting waves with substantially smaller amplitudes. 
But in any case, coupling coefficients are of order /c 7 with 7 > which apparently con- 
tradicts to the existence of Kolmogorov-like power spectra k' 7 with 7 < in the region of 
short waves. 



4. Laminated WT 

In order to resolve this apparent contradiction we have go back to the very base of the 
KAM theory. Its main results are based on the famous Thue theorem, giving low estimate 
for the distance between any algebraic number a of degree n > 2 and a rational number 
p/q e Q: 

P, c(a) 
\a--\> A _V _/ /9 , Ve>0 

where c(a) is a constant depending on a and A can be arbitrary small. This fact allows to 
construct KAM tori, with a being a ratio of frequencies of interacting waves, and KAM 
theorem states then that almost all tori are preserved. "Almost" means in particular 
that tori with rationally related frequencies (corresponds to a being an algebraic number 
of degree n = 1) are explicitly excluded from consideration. Since the union of invariant 
tori has positive Liouville measure and Q has measure 0, this exclusion is supposed to be 
not very important. 

Coming back to the example of spherical planetary waves, one can see immediately 
that the ratio of the frequencies for resonantly interacting waves in this case is a rational 
number, both for exact and approximate interactions. Therefore, these waves are not 
described by KAM theory. Obviously, the same keeps true for an arbitrary wave system 
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with rational dispersion function uj. Case of planetary waves in a square domain is a bit 
more complicated. It is proven jllj that exact resonances in this case are described by 
the wave vectors ki = (mj,nj) satisfying the (necessary) condition 

h = a iy /q, at £N Vi = 1, 2, 3, 

with the same square-free q. This fact allows to construct disjoint classes of resonantly 
interacting waves and q is called index of the class. Obviously, for the waves belonging to 
the same class, the ratio of their frequencies is a rational number, uJi/coj £ Q, and 
these waves are excluded from KAM theory. Waves, interacting approximately, may have 
different indices (not necessarily) but the ratio of frequencies uJi/ujj is then an algebraic 
number of degree < 2 due to the form of dispersion function and, therefore, these waves are 
also excluded from the KAM theory. Similar results can be proven for exact resonances 
in many wave systems, for instance, for an arbitrary wave system in which dispersion 
function is a polynomial of finite degree on k with at least one non-zero coefficient in front 
of an odd degree of k. 

Let us summarize the results obtained. Continuous WT (CWT) describes energetic 
behavior of a wave system for the short-waves' part of spectrum excluding nodes of 
rational lattice thus leaving some gaps in the spectrum which are supposed to be not 
important in short-waves' part. Discrete WT (DWT) fills these gaps all over the spectrum. 
In fact we have two layers of turbulence - CWT (layer I) and DWT (layer II), which are 
mutually complementary and should be regarded simultaneously. 

Layer I provides KAM tori and stochastic enough turbulence in short-waves range 
with Kolmogorov's spectra in the inertial interval; direct or inverse energy cascades are 
possible; wave-numbers range of energy pumping influences the results. 

Layer II provides a countable number of waves with big amplitudes all over the wave 
spectrum; some of the waves do not change their energies (non-interacting waves) and 
others do exchange energy within small independent groups; there is no energy cascade 
at this layer; results do not depend on the wave-numbers range of energy pumping. 

The co-existence of these two layers means in particularly that in the waves' range, 
classically described by Kolmogorov-like spectra, there exist also a small group of waves 
from layer II with substantially bigger energies than their "neighbors" from layer I, and 
this small group generates appearance of some structures. We give here a few examples of 
known phenomena which can be explained in the same frame of the model of laminated 
turbulence. 

(1) Very clear example of the co-existence of these two layers is given in ^Jj where tur- 
bulence of capillary waves was studied in the frame of simplified dynamical equations for 
the potential flow of an ideal incompressible fluid. A stationary regime of so-called "frozen 
turbulence" had been discovered: in small wave-numbers region wave spectrum consists 
of "several dozens of excited low- number harmonics" which construct "ring structures in 
the spectrum of surface elevation" . The appearance of these structures does not depend 
on the damping and pumping, and in all computations "the Kolmogorov's spectrum co- 
exists with the spectrum of another, "frozen" type, concentrated in the region of low 
wave-numbers and fastly decreasing to large wave-numbers. If the level of nonlinearity is 
low enough, such "frozen" regimes are dominant" (|17j). Obviously, these ring structures 
are due to non-interacting waves of layer II and similar structures were also observed in 
laboratory experiments and identified as such Notice that as there exists no exact 
three- wave interactions among capillary waves with u> 2 = k 3 we observe in regime of 
frozen turbulence of capillary waves only discrete waves with constant amplitudes. 

(2) Similar experiments/simulations with, say, four- waves interactions among gravity 
waves, to 2 = k, will demonstrate that frozen turbulence partly "thaws out" because also 
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changes in the amplitudes of resonantly interacting discrete waves should be observed (cf. 
"bursty" spectrum in [T5]). 

(3) Mesoscopic turbulence (corresponds to "transition" interval mentioned above) 
discovered in numerical experiments on modelling of turbulence of gravity waves on the 
surface of deep ideal incompressible fluid gives another example of manifestation of lam- 
inated turbulence. It was established that not all waves "have the same rights" and the 
existence of so-called "elite society of harmonics" has been demonstrated; their number 
amounts to only 6% of the total number of harmonics being 10 4 but they play the most 
active role in mesoscopic turbulence. These elite harmonics correspond to exact and ap- 
proximate resonances of the discrete waves (layer II). The number of these harmonics 
depends, of course, on the specific wave system. For instance computations with spherical 
planetary waves in the domain of wave numbers < m, n < 1000 (which is far beyond the 
region of applicability of BVE) shows that total number of harmonics taking part in exact 
resonant interactions is 22683, i.e. slightly more than 2%. Consideration of approximate 
interactions of the layer II increases this amount but not substantially. 

(4) Zonal extended vertices (flows in latitudinal direction) in the atmosphere can 
possibly be explained in terms of plane images of spherical waves with coupling coefficient 
n 3 / 2 . 

At the end of this letter we would like to make one important remark. To describe the 
short waves of layer II it is necessary to develop fast algorithms of solving Diophantine 
equations in very big integers (of the order 10 12 and more). We consider it possible basing 
on the existence of disjoint classes of waves participating in a single solution. This is our 
current object of interest. 

Author is sincerely grateful to V.E. Zakharov for his encouragement and permanent 
interest for this work, and whose numerous results in the theory of wave turbulence gave 
us necessary insights while creating a model of laminated turbulence. Author is also very 
much obliged to anonymous referees for their valuable remarks and suggestions. 
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